Measuring the momentum of a nanomechanical oscillator through the use of two 

tunnel junctions 
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We propose a way to measure the momentum p of a nanomechanical oscillator. The p-detector is 
based on two tunnel junctions in an Aharonov-Bohm-type setup. One of the tunneling amplitudes 
depends on the motion of the oscillator, the other one not. Although the coupling between the 
detector and the oscillator is assumed to be linear in the position x of the oscillator, it turns out 
that the finite-frequency noise output of the detector will in general contain a term proportional 
to the momentum spectrum of the oscillator. This is a true quantum phenomenon, which can be 
realized in practice if the phase of the tunneling amplitude of the detector is tuned by the Aharonov- 
Bohm flux "1> to a p-sensitive value. 

PACS numbers: 72.70.+m,73.23.-b,85.85.+j 
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In nano-electromechanical (NEM) systems, it is the po- 
sition of the oscillator that typical measurement devices 
(like tunnel junctions or single electron transistors) are 
coupled to. Using these detectors, position measure- 
ments with sensitivities close to the standard quantum 
limit have already been observed [H, [2j S| • From a funda- 
mental point of view, it is desirable to go further, i.e. to 
prepare and manipulate NEM oscillators in the quantum 
regime. A quantum NEM system would allow us to study 
an ideal realization of a continuous variable quantum sys- 
tem The exploration of such systems has to be seen 
as complementary to the wide study of two-level systems 
done in the context of quantum computing. 

In order to be able to fully characterize a continuous 
variable quantum system that is described by two non- 
commuting operators x and p, we need to be able to mea- 
sure expectation values of moments of both of them Q . 
Only this allows, for instance, to detect the entanglement 
between two (or more) NEM devices J^. The literature 
already contains proposals regarding quantum measure- 
ments of the momentum of macroscopic objects like those 
used for gravity- wave detection Q. However, none of 
these proposals have been realized in practice. In this Let- 
ter, we propose a realistic and feasible way to measure the 
momentum of a nanometer- sized resonator. This is a non- 
trivial task since the coupling between the detector and 
the oscillator is naturally described by an x-dependence 
but not a p-dependence. Nevertheless, the proposed 
setup (shown in Fig. [T]d) allows for a measurement of 
the momentum spectrum Sp{u}) — J c?te"^'({p(t),p(0)}) 
of the oscillator. This can be done because we have 
found a way to tune the phase of the tunnel coupling 
term that is sensitive to the position of the oscillator 
by an Aharonov-Bohm (AB) flux see Fig. [T}3. Re- 
lated setups have been investigated recently in the con- 
text of dephasing due to the coupling of an AB ring 
structure to a NEM device Q . A typical position detec- 



tor which has been analyzed theoretically in great detail 



^10,11, [13, M , [Hi and experimentally realized 
la . 17 1 is depicted in Fig. [T^. It shows a single tun- 
nel junction coupled to a NEM oscillator. A thorough 
analysis of the coupled quantum system leads to the re- 
sult that the output signal of the detector is sensitive to 
the position spectrum Sx{o-!) = J dte"^'({a;(i), i;(0)}) of 
the oscillator. The modification of the detector shown in 
Fig. [lb instead allows for a measurement of Sp{Lu). 
The Hamiltonian of the coupled system H = i?osc + Hb + 
Htun is the sum of the Hamiltonian of the (quantum) 
harmonic oscillator i?osc (with mass M and frequency f2), 
the bath Hamiltonian Hb (describing the leads of the 
detector), and the tunneling Hamiltonian i?tun (which 
couples the dynamics of the electrons that tunnel across 
the junction to the motion of the oscillator): 
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Here, k (q) is a wave-vector in the right (left) lead, c*^^^ is 
the electron annihilation (creation) operator, and F^t) 
is an operator that decreases (increases) m, the num- 
ber of electrons that have tunneled through the system, 
by one. It allows one to keep track of the transport 
processes during the evolution of the system. For typ- 
ical nanoresonators, the mass M varies between 10~^^ 
and 10~^^ kg, and the resonance frequency is usually 
IMHz <n/2Tr < IGHz 

We first discuss the model in the standard configuration 
shown in Fig. [iji and later on describe the new setup 
in Fig. [Hd. For small displacements with respect to the 
tunneling length (which is the relevant regime in typi- 



FIG. 1; (color online) a Position detector. The figure shows schematically a position detector of the motion of a NEM oscillator 
(red). The detector is based on a tunnel junction with a tunnel matrix element T{x) which depends on the position of the 
oscillator. The shaded regions (yellow) are assumed to be conducting, b Position and/or momentum detector. The figure 
illustrates a detector which contains two tunnel junctions that form a loop threaded by a magnetic fiux <E>. Tuning the fiux will 
change the performance of the detector from being able to detect the power spectrum of the position operator of the oscillator 
to being able to detect the power spectrum of the momentum operator of the oscillator. For clarity, the two insets show a 
simplified illustration of the detectors in a and b. 



cal experiments on NEM devices), the tunneling ampli- 
tude T{x) can be taken as a linear function of x, namely 
r(^) = (e*'^" /27rA) (tq + e^^'ni) , where tq and n are 
real, and A is the density of states. The phases ipo and 
T] describe details of the detector-oscillator coupling [l^. 
It can be shown that our detector is quantum-limited for 
any 77, according to the definition of a quantum-limited 
detector in Ref. Within the single-junction setup 



(Fig. [T^), the relative phase r] is sample-dependent and 
cannot be tuned experimentally. In a typical device, 
the x-dependence of the phase of T{x) is much weaker 
than the cc-dependence of the amplitude. Then, we can 
set ?7 ~ 0, and the tunnel junction acts as a position- 
to-current amplifier where the frequency-dependent cur- 
rent noise Si{lo) of the detector contains a term propor- 
tional to the position spectrum Sx{uj) of the oscillator, 
i.e. ASi{uj) = St(lo) — 2e{I) k \l.Sx{oo) with the gain 
of the amplifier [isl |20| . 

We now demonstrate that a tunnel junction with a 
phase 77 — 7r/2mod7r acts as a momentum detector 
and ASi{lj) w XpSp{uj), where Xp is the gain of the 
momentum-to-current linear amplifier. The critical re- 
quirement to build a momentum detector is to be able 
to vary rj experimentally. This can be done using the 
AB-type setup shown in Fig. [T)d: a metallic ring where 
one arm is a standard tunnel junction position detector 
with tunneling amplitude Td{x), and the other arm is a 
position- independent tunnel junction with tunneling am- 
plitude T„ 2l|. The total transmission amplitude T(x, $) 
of the device is the sum of both tunneling amplitudes [23| . 
Since only one arm shows a position-dependence, the in- 
duced phase difference between the two arms affects the 
position-independent and the position-dependent parts 
of the tunneling amplitudes tq and ti in a different way. 
Explicit calculation shows that the tunneling amplitude 
is given (up to a global gauge-dependent phase factor) by 



r(x,$) = (to($) -^-e'''(*)Ti£)/(27rA) with 



^o(*) = '^o.d + '^0% + 2T-o,dTo,„cos(^27r— + (/3o,d - <i2o,«), 



Vi"^) = 27r— + (pi,d - lpq. 
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where we have defined T, 



ro,„, Td EE e"^"-ro,d ■ 
e"^'--''Ti^dX, Ti = Ti^d, and $0 = h/e 24]. The position- 
independent part of the tunneling amplitude to($) dis- 
plays the standard AB oscillations as a function of flux. 
Likewise, the relative phase ??($) shows a distinct de- 
pendence on the flux. Importantly, for tq^u > 'Ha.d, 
the phase ??($) can be tuned continuously in the whole 
range [— 7r,7r]. In the limit, where ro,„ 3> TQ di vi^) ^ 
27r^ -t- ri{^ = 0) varies linearly with the applied flux. 
In the opposite regime tq.u ^ To.d, rj no longer depends 
on <I>. Therefore, it is crucial to put the tunneling am- 
plitudes in the regime where ?7($) can be tuned to n/2. 
We will show below that a feasible way to calibrate ri{^) 
to the p-sensitive point tt/2 is a measurement of the flux 
dependence of the current through the AB detector. 
Similarly to Refs. l^, 20j], we study the coupled system 



using the quantum equation of motion for the charge- 
resolved density matrix within the Born-Markov approx- 
imation, assuming that eV 3> hfl. 

It has been derived previously that, under the assump- 
tion that the tunneling amplitude depends linearly on x, 
the equation of motion for the reduced density matrix 
of the oscillator is of Caldeira-Leggett form [l^j [Si ■ 
Thus, it contains both a damping and a diffusion term. 
When the electron temperature is much smaller than the 
applied bias (and taking V > 0), the detector-induced 
damping coefficient is 7+ = Ht^ / (iirM) and the diffusion 
coefficient is D+ = 2Mj+kBTcs with Tcs = eV/2kB. 



3 



In general, the oscillator is not only coupled to the de- 
tector but also to the environment. The coupling to this 
additional bath is controlled via 70 = ^/Qo (related to 
the finite quality factor Qo of the mode, which in current 
experiments varies from 10"^ to 10^ [18|) and the associ- 
ated diffusion constant = 2MjokBTenv that must be 
added to the detector-induced damping and diffusion con- 
stants to find the total damping coefficient 7tot = 7+ +7o 
and the total diffusion coefficient Dtot — Dq + . Tcnv 
denotes the temperature of the environment. In typical 
experiments, it varies from SOmK to lOK. Within our 
model, all the these system parameters are independent 
of the applied flux. 

It is now straightforward to calculate the current and 
the current noise of the detector. We skip the details 
here (see Ref. for 77 = 0) and directly turn to the 
results. The average current of the detector is given by 



trum is therefore AS* = ASi + AS2 with 



I = — {t§ +2cosi]ToTi{x) + tI{x'^)) 



h 

2e7+ro 



sin 77 (p) — 67+. 
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For rj ^ mod n, the average current contains a term pro- 
portional to the average momentum of the oscillator that 
does not vary with the applied bias fl^. However, since 
(p) = in the steady-state, the average current contains 
no information about the momentum of the oscillator. 
Therefore, the current of the detector can not be used as 
a p-detector in the steady-state. Nevertheless, the cur- 
rent is important to calibrate 77 to the p-sensitive value 
tt/2. a careful analysis of the current / as a function of $ 
shows that the inflection points of /($) correspond pre- 
cisely to values of 77 = 7r/2mod7r. Therefore, we can use 
a current measurement to tune 7/ to a p-sensitive value. 

In the experimentally relevant regime, where T^(a;^) ^ 
Tq, and for uj ft, the dominant contributions to the 
current power spectrum of the detector are 



Si{uj) = 2e{I) + 8e^cj / dt sin(cjt) x 
Jo 

^ cos Tyro Ti((a;77i)) - sin7;((p777)) 



(6) 



where {{ab)) = {ab) — {a){b). We now further analyze 
the added noise due to the presence of the oscillator, 
AS ~ Si{uj) — 2e{I). This noise spectrum is the sum 
of a contribution arising due to correlations between the 
transfered charge 777, and position (term ~ {{xm}) in Eq. 
([6])), which we call AS*!, and one due to correlations 
between 771 and the momentum of the oscillator (term 
~ {{pm)) in Eq. which we call AS'2. The full spec- 
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where the position and the momentum gain are given 
by Xx = 2eToTi{eV/h) COST/ and Xp = (e/27rM)Tori sin 7;, 
respectively. We now discuss several limits of the cur- 
rent noise Si{uj) of the detector in the case of a gen- 
eral phase rj. For 77 — OmodTr, we recover Eq. (30) of 
Ref. [2^ - the position detector result. More interest- 
ingly, for 77 — 7r/2mod7r, A^; = and the detector output 
contains only two terms: The flrst one is proportional to 
Sp{uj) and therefore peaked around il. The second one is 
proportional to (fJ^ — cj^) and contributes negligibly near 
resonance u fl. Hence, for 7/ = 7r/2, we obtain 



ASiu; « f]) « A" 1 - 



MeV \ 
W))) 



Spiu). 



(9) 



Thus, the added noise is directly proportional to the mo- 
mentum spectrum of the oscillator. This is the key result 
of our Letter. 

From the parameter dependence of each gain, we can es- 
timate that the momentum signal at 77 = 7r/2 should be 
typically smaller than the position signal at 77 = by a 
factor {eV/hV,)^. Nevertheless, it is unambiguously possi- 
ble to identify ap signal in the current noise. We now de- 
scribe three different ways to do this. First, since Xx (xV 
while Ap is independent ofV, the bias voltage dependence 
of the noise spectrum can also be used to confirm that 
momentum fluctuations are measured. Secondly, for an 
oscillator undergoing Brownian motion, the temperature 
dependence of both signals differs qualitatively. Like in 
the position detector case, the momentum signal is re- 
duced by a quantum correction (the term proportional to 
—MeV/dp^)) in Eq. ©) that arises from the finite com- 
mutator of X and p. However, there is a fundamental dif- 
ference between the x-detector result (Eq. (7) of Ref. [ist ) 
and the p-detector result (Eq. ([9])). In the former case, 
the quantum corrections are always small compared to 
the leading terms and therefore the peak at resonance 
is always positive. In contrast, the two terms in Eq. ([9]) 
can be of equal magnitude and compete about the sign 
of AS{uj w il). The p-sensitive current noise in Eq. ^ 
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changes sign when the effective temperature of the oscil- 
lator is equal to (eT^/fc_B)(l-7+/27tot)/(l-7+/7tot)- For 
a cold environment Tonv ^ eV^, AS{uj) is negative at the 
resonance, whereas, for a hot environment T^^v > sV, 
AS{u!) is positive. This change of sign never appears 
during a position measurement, so this pronounced dif- 
ference between a x-dependent and a p-dependent signal 
can be used distinguish the two. We illustrate the change 
of sign in the inset of Fig. [21 where the added current 
noise for 77 = 7r/2 is plotted for different Tonvll^- 



We have demonstrated how the proposed detector can 
be made sensitive to either displacement or momentum 
fluctuations of the oscillator. 
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FIG. 2: (color online) Added current noise (normalized by 
2e/o = 2e^ToV/h) of the proposed momentum detector due to 
the presence of the oscillator. For all curves, the bias is eV = 
SOftn, 7+ — 7tot/4 and 7tot = $1/200. The main panel shows 
the total detector output for different values of the tunneling 
phase rj and for Tcnv = 0. The (blue) solid, (black) dotted, and 
(red) dashed lines correspond to 2r]/n — 1, 1.005, and 1.01, 
respectively. In the inset, the current noise at the p-sensitive 
phase 77 — n/2 is plotted for two different temperatures of 
the environment Tenv ~ (solid line) and Tcnv = beV/ks 
(dash-dotted line). 

In the main panel of Fig. [21 we plot the full detector out- 
put for different values of i] near the optimal operation 
point for momentum detection. Away from 77 = 7r/2, con- 
tributions to the current noise become important 
and wash out the momentum signal ~ A^. Indeed, for 
small A77 = rj — 7r/2, the ratio \x/\ ^ {—lS.r]){eV/Ml) 
of the two amplification factors becomes large as soon as 
jAryj > Ml/eV. In the high-bias regime {eV ^ fifi), mo- 
mentum detection therefore requires good experimental 
control over the applied flux. At moderate bias eV > Ml, 
the requirement on A?7 becomes less restrictive. Finally, 
the current noise spectrum at rj = 7r/2 shows a strong 
symmetry around that makes the optimal operation 
point easily identifiable. 

In conclusion, we have shown how a modified tunnel junc- 
tion position detector can be designed to detect the mo- 
mentum fluctuations of a NEM oscillator. By using two 
tunnel junctions in an AB-type setup, it is possible to 
precisely tailor the interaction Hamiltonian between the 
detector and the oscillator via an external magnetic field. 
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